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Abstract 

We investigate the solar-system constraint on the f{R) theory of modified gravity with chameleon 
mechanism, where f(R) represents the deviation from general relativity in the gravity action. We 
obtain a stringent bound to a general, non-constant deviation function f(R). — 10~ 15 < df/dR < 
when R/Hl ~ 3 x 10 5 , and a loose bound: < Rd 2 f/dR 2 < 2/5 when R/H$ > 3 x 10 5 , by 
requiring the thin-shell condition in the solar system, particularly in the atmosphere of the Earth. 
These bounds can be conveniently utilized to test the f(R) models with given functional forms of 
f(R) and to obtain the constraints on the parameters therein. For demonstration we apply these 
bounds to several widely considered f{R) models. 
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I. INTRODUCTION 

The discovery of the cosmic acceleration at the present epoch [l,|2| indicates the existence 
of repulsive gravity that dominates the present universe at large scales. It may be an 
indication of an energy source of anti-gravity or a sign of the modification of the gravity 
theory. Along the latter direction the f(R) theory of modified gravity (for a review, 



see [14|) has been proposed as a possible explanation of the cosmic acceleration, where the 
gravity is described by a function of the Ricci scalar R in the action. On the other hand, 
the f(R) theory can be treated as a simple way of modeling the possible deviation from 
Einstein's general relativity (GR). 



As an essence of cosmology, f(R) gravity needs to pass the cosmological test [15|-|24| 
involving the observations of the cosmic expansion and the cosmic structure formation. 
As a gravity theory, it needs to pass the local gravity test that in general gives the most 



stringent constraint on modified gravity so far 



m 



16 



25 



4281] . It has been pointed out 



that f(R) gravity, as a special case of the scalar-tensor theory, may pass the solar-system 
test with the help of the "chameleon mechanism" 0, [k]] , meanwhile driving the late-time 
cosmic acceleration. In the chameleon mechanism the scalar field of the scalar-tensor theory 
can behave differently in different environments, depending on the ambient mass density. 
This feature makes it possible to have significant deviations from GR at the cosmological 
scales at late times and meanwhile have tiny deviations both in the solar system at present 
and at all scales at early times. 

In this paper we investigate the solar-system constraint on the f(R) gravity with the 
chameleon mechanism. 1 As a result, we obtain the following constraints on a general, non- 
constant function f(R) that represents the deviation from GR in the gravity action. 

- 10~ 15 < df/dR < when R/H% ~3x 10 5 , (1) 
< Rd 2 f/dR 2 < 2/5 when R/H 2 > 3 x 10 5 . (2) 

These constraints will be derived in Sec. IIVI Before that, in Sec. HH we will introduce f(R) 
gravity, the chameleon mechanism and the thin-shell condition, and in Sec. II III elaborate on 
the thin-shell parameter in f(R) gravity and obtain its relation to f{R), with which the 
solar-system constraint on the thin-shell parameter can be transferred to that on f(R). 



1 We consider the metric formalism of f(R) gravity. 
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II. f(R) GRAVITY WITH CHAMELEON MECHANISM 



We consider the f(R) theory of modified gravity with the action, 2 
S = 2^ / dAx ^9 [R + f{R)) + S m ( 9fMV , * m ) , 



(3) 



where f(R) is a function of the Ricci scalar R and represents the deviation from GR, k is the 
gravitational constant, S m the matter action and \l/ m the matter field. This theory can be 
transformed to a scalar-tensor theory in the Einstein frame via a conformal transformation, 



9iw 



:i + f R )g,u = e- 2 ^ gflu , /3 = -l/v / 6 ; 



(4) 



where f R = df/dR. The resultant action in the Einstein frame is 

1 ^ 1 



S f 



d x\J — g 



S m (e 2 ^% v ,V m ) 



(5) 



where 



1 



\n[l + f R (R)} 



V(<f>) 



2/3k 

Rf R (R) - f(R) 



2k* [1 + f R (R)] 2 ' 

R = R[g»u\. 



(6) 

(7) 
(8) 



The notation about the derivatives of f(R), Jr = df/dR and f RR = d 2 f/dR 2 , will be used 
in the remainder of this paper. 

The potential V in Eq. (j7]) is a function of the scalar field 0, provided that the Ricci 
scalar R is a function of 4> given by Eq. (E]) with a well-defined inverse function f^ 1 . In 
the action Se the scalar field acquires an additional coupling to matter through the Jordan- 
frame metric tensor g^ v in the matter action S m . This additional coupling is the key to the 
chameleon mechanism. It makes the scalar field <fi behave differently in different environ- 
ments, depending on the ambient mass density. 

To study the f(R) chameleon gravity with the action Se, for simplicity we consider 
a spherically symmetric system, especially a sphere with different constant mass densities 
inside and outside the sphere. In the following we will present the field equation and the 



We use the convention, {- 



-}, for the metric signature. 
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solution for the scalar field <f>(r), where r is the physical distance from the center of the 
system. 

For a spherically symmetric space-time with a mass density distribution p*(r), the field 
equation is 

d 2 <f> | 2d<t> _ dv eS (4>, P *) 

dr 2 r dr d(f) 

where the effective potential 

V eS (cf>,p*) = V(<f>) + e^p*. (fO) 

The mass density p* is a conserved quantity in the Einstein frame. By p* = e 3l3K ^p, it is 
related to the mass density p in the f(R) gravity with the action S. 

To satisfy the constraints from the solar-system experiments, the effective potential V e g 
should have a minimum and should be steep around the minimum when the mass density p* 
is equal to those in the solar system, including the regions inside the Earth (p^ ~ 5.5 g/cm 3 ), 
in the atmosphere (p* tm ~ 1CT 3 g/cm 3 ) and outside the atmosphere (p G ~ 10~ 24 g/cm 3 ). 
Let (p m denote the location of the potential minimum and m the mass of 0, i.e., 

fyKff(0m,P*) =0, (11) 

m 2 (p*) = <9 2 V^0 m ,p*). (12) 

The basic condition raised above is then recast by the requirements: (1) the existence of 
m , (2) positive m 2 , and (3) large m (i.e. short Compton wavelength), when p* ~ p^, p* tm , 
p* G , etc. 

We consider a sphere with the radius r s and with different constant mass densities, p* n 
and p* ut , inside and outside the sphere. The locations of the potential minima C and 0^, 
the masses m in and m out , the total mass of the sphere M s and the Newtonian potential <3> s 
at r s (i.e. on the surface of the sphere) are defined as follows. 

d</>V eS (0 C , p* n ) = 0, dtVrf (0oo, p* out ) = 0, (13) 

< = dlV cS (0 C , p* n ) , m 2 ut = dlV cS (0^, plJ , (14) 

„ . 4:71 c, K 2 M s 

M, S -rlpl, ». SsF -. (15) 
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The field profile, solution of Eq. ([9]), is 

cf) (r < r a ) ~ 4> c , 



(r>rs) " ~ w; v^T; r ' 

when the following thin-shell condition is satisfied. 

< ^ s ^2 - 6j « i. (17) 

The thin-shell parameter Ar s /r s is proportional to the ratio of two potential differences, 
k^oo — K(p c and $ s — [where $oo = $(r = oo) = 0], which respectively relate to the 
strength of the fifth force induced by <fi and that of the Newtonian gravitational force. Thus, 
roughly speaking, the thin-shell condition requires the weakness of the fifth force compared 
to the Newtonian gravity. 

The profile in Eq. (fT6|) depends on the thin-shell parameter Ar s /r s that is determined 
when the function f(R) and the mass densities, p* n and p* ut , are given. Accordingly, the 
solar-system bounds to the fifth force can constrain Ar s /r s and thereby constrain f(R). 

III. THIN-SHELL PARAMETER IN f(R) GRAVITY 

Here we will derive the relation between the thin-shell parameter Ar s /r s and the function 
f(R)- With this relation at hand we can obtain the constraint on general f(R) from the 
solar-system tests of gravity which give bounds to the thin-shell parameter. In addition, the 
requirement of the existence of the potential minima, mf a > and ml nt > 0, also gives a 
basic condition of f(R), for which we will derive the relation between f(R) and the masses, 
m in and m out . 

The thin-shell parameter is related to f(R) through the locations of the potential minima, 
C and 0oo, which satisfy Eq. QEJ). With dV/d<p = (dV / dR)(dR/ defy and = l/y/l + f R 
we have 

«9 K ff (0,p*) = ^L + fiKePxtp* (18) 
R + 2f-Rf R | p K p* 
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d 2 V 
d(j> 2 
1 



3R + 4f-Rf R + I3 2 k 2 p * 



(20) 
(21) 



' 3f RR 3(1 + f R ) 2 ' VTTh 
In the solar system the deviation from GR must be small, i.e., \ f\<^R and \f R \ <C 1, which 



we will use for the approximation involved in the following derivation. Therefore, 

1 



[R ~ «V) + (2/ - 3Rf R + ^/r«V) 



dlV cS 



-±--r+ \k 2 p *) + (~f + 7 -Rf R - y R K 2 P * 

ijRR 6 / \ 3 3 12 



(22) 



(23) 



The location of the potential minimum, m , is given by c^T4ff(0 m , p*)—0, and therefore 
satisfies 

[R + 2f-5Rf R /2L~ K y, 



or, to the lowest order, 



Thus, 



R{<f> m )~K 2 P * 



<f*r, 



In [1 + f R (R ~ K 2 p*)] 



2(3k 
1 

'2/3k 



f R {R~K 2 p*). 



We then obtain the formula for the mass, 
m 2 (p*)=^y efr (0 m ,p*)~ 



and the relation between the thin-shell parameter and f(R), 



(24) 
(25) 

(26) 
(27) 

(28) 



Ar, 



Pout) - /r(« P*n) 



r s " 2$ s 



(29) 



In many viable /(-R) models, |/j?(i?i)| 3> \f R {R2)\ when i?i <C i? 2 , in order to fit the 
cosmic microwave background (CMB) and the big-bang nucleosynthesis (BBN) observational 
results that require the deviation from GR be tiny at early times, meanwhile generating the 
cosmic acceleration with a significant deviation from GR at late times. In this case, 



Ar s /r s ~ -f R (K 2 p* out )/2<f> s when p* out < p* 



(30) 



The condition p* ut <C p* n is satisfied in many solar-system experiments. 



IV. SOLAR-SYSTEM CONSTRAINTS ON f(R) GRAVITY 



To fit the solar-system constraints, the thin-shell condition in Eq. (iT7|) needs to be satisfied 
in the solar system. This condition requires (1) the existence of the minimum of the effective 
potential V e s(cj),p*) and (2) the smallness of \ fn\. 

From Eq. (1281) . the existence of the minimum of V e s, i.e. m 2 > 0, entails (to the lowest 
order) 

< Rf RR < 2/5 when R ~ k 2 p* > 3 x 10 5 # 2 , (31) 

where we have considered various environments in the solar system, with the mass density p* 
ranging from 10~ 24 g/cm 3 (Space) to 5.5 g/cm 3 (Earth) and accordingly with R/Hq ranging 
from 3 x 10 5 to O(10 30 ). For simplicity we use R > 3 x 10 5 Hq in the above expression. 

Regarding the thin-shell condition, an upper bound (Ar s /r s ) max of the thin-shell param- 
eter constrains f R via Eqs. ffT7|) and fl29|) as follows. 

< f R (K 2 pl) - f R (K 2 p* ont ) < 2$ s • (Ar s /r s ) max « $ s , (32) 

which gives an upper bound to the f R variation. In the cases where 

|/i?(« 2 Pm)| < l/flfaVout)! f O r Pm > Pout , (33) 



it gives an upper bound to |/r(k 2 p* 



out/ 



< -f R < 2$ s • (Ar s /r s ) max « $ s when R~ K 2 p* ont . (34) 

In the following we will consider the upper bound (Ar s /r s ) max obtained from the solar-system 
experiments involving the Sun and the Earth, respectively. 



The experimental tests of the post Newtonian parameters in the solar system 



27] 



3l| give 



(Ar s /r s ) max = 1.15 x 1(T 5 . (35) 
With $ ~ 2.12 x 10~ 6 for the Sun and p* ut ~ p* G ~ 10~ 24 g/cm 3 , we obtain 

- 5 x 1CT 11 < f R < when R/H 2 ~ 3 x 10 5 . (36) 

The case of the Earth experiments is more complicated. In this case the spherically 
symmetric system has three regions: the Earth, the atmosphere, and the beyond, with the 
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p*{r) = < 



mass densities as follows. 

~ 5.5 g/cm 3 for < r < r e , 

Patm - 10_3 g/ Cm3 fOT r ® < r < r atm , ( 37 ) 

p* G ~ 10" 24 g/cm 3 for r > r atm , 

where the Earth radius ~ 6.4 x 10 3 km and the thickness of the atmosphere Ad atm = 
atm — — 10-100 km. With the following thin-shell condition satisfied, 
Ar atm _ K (<f> G - atm ) 1 , 2 

J- , ^atm = - K Patm r atm 5 l d °J 



'"atm 



the scalar field profile is [301 ] : 



6/3$atm ' 6 



</>(r) 



''atm 



(3k\ /3Ar e \ M e e- mG ( r - ratm ) 

h G for r > r atm 



for < r < r e , 
for r® < r < r atm , (39) 



47T / \ r g 

Ar e _ - atm ) _ 1 ^ _ 1 . o ,. .2 



where 0®, atm and 0g respectively denote the locations of the effective potential V e g minima 
in the three regions, M® the mass of the Earth, and and $ a t m the Newtonian potentials: 
M e ~ 6 x 10 24 kg, $ e ~ 7 x 10~ 10 and $ atm ~ 10~ 13 . Note that the shin-shell condition 
in Eq. fl38|) automatically leads to Ar ffi /r® <C 1. The experimental bounds of the thin-shell 
parameters, < Ar»/ri < (Arj/rj) max <C 1 for % = ©,atm, then give upper bounds to |/#|: 

< -f R (K 2 p G ) < 2$ ?; • (An/r^^ < , 2 = ©, atm , (41) 

in the case where |/ R («: 2 p atm )| < \ f R {^ 2 p* G )\ for p* atm > 

An essential experimental bound comes from the basic requirement that the atmosphere 
has a thin shell, i.e., the thickness of the thin shell should be smaller than that of the atmo- 
sphere: < Ar atm /r atm < Ad atm /r atm . Taking the thickness of the atmosphere Ad atm — 50 
km, i.e. Ac/ atm /r atm ~ 8 x 10~ 3 , we obtain a very stringent bound to 

- 10" 15 < f R < when R/Hl ~ 3 x 10 5 . (42) 

To sum up, from the thin-shell condition in the solar system we have obtained a constraint 
on Rfun in Eq. (I3T]) and a stringent upper bound to \f R \ in Eq. (142]) for non-constant f(R)- 3 



3 Note that fa and fan can vanish in the case of constant f(R) that is equivalent to the ACDM model and 
therefore not taken into consideration when we investigate f(R) gravity for possible deviations from GR. 
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They respectively come from the requirement of the existence of the effective potential 
V e ff((j), p*) minimum in the solar system and from the thin-shell condition in the atmosphere, 
with the precondition that |/| <C R and \fn\ <C 1 in the solar system, and |/R(^ 2 p* tm )| <C 

\f R (K 2 ph)\- 

The constraints on f(R) we obtained can be conveniently applied to the f(R) models 
where the functional forms of f(R) are given. For demonstration, here we apply the con- 
straint on fn in Eq. (|42p to the following widely considered models 20]: 



f(R) = -\RJ 1 (x), x = R/R c , X,R c >0; (43) 

(1) fi(x) = xP, 0<p< 1; 

(2) Hu and Sawicki Q: fi{x) = x 2n /{x 2n + 1), n > ; 

(3) Starobinsky fl: f x (x) = 1 - (1 + x 2 )"", n > ; 

(4) Tsujikawa [h]]: fi(x) = tanh(a;) ; 

(5) Linder Q: fi(x) = 1 - . 

Considering A ~ -Rc ~ ^ 2 Pc (where p c is the critical density at the present time) 

and R/R c ~ 10 5 , we find that the constraint in Eq. (l4"2l) requires p < 10~ 10 in Model 
(1) and n > 1 in Models (2) and (3), and it is well satisfied in Models (4) and (5). In 
addition, the constraint on RJrr in Eq. (!3T!) is also satisfied in these models. Note that 
in these models and under the above consideration, the conditions |/| <C R, \fn\ 1 and 
|//j(K 2 p* tm )| |/r( k2 Pg)I are satisfied, so that we can legitimately use the simple and 
stringent constraint on /^(k 2 Pq), i.e. the constraint in Eq. ( )42l) . 



V. CONCLUSION 



The cosmological tests and the local tests give essential constraints on the deviation 
from GR in f(R) gravity at different values of the Ricci scalar R, which correspond to the 
constraints at the cosmological scales at different epochs of the cosmic expansion history. 
Roughly speaking, CMB and BBN stringently constrain the large-i? (early-time) behavior 
of f(R), the solar-system tests stringently constrain the moderate- R (middle-age) behavior, 
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and the cosmological observations about the late-time universe constrain the small- R (late- 
time) behavior. 

With regard of the cosmic history, the constraints on the f(R) modified gravity at the 
cosmological scales from the early times to the present are summarized as follows. 

• z > 10 3 : The deviation from GR, such as f R , should be small at early times when the 

redshift z > 10 3 , as required by the CMB and the BBN observations. 

• z > 70: < R/rr < 2/5 when R/Hq > 3 x 10 5 that roughly corresponds to the epoch 

z > 70. This is a basic requirement in the chameleon mechanism for the solar system. 

• z ~ 70: — 10~ 15 < fjt < when R/Hq ~ 3x 10 5 that roughly corresponds to the time when 

z ~ 70. We obtain this stringent constraint from the thin-shell condition required in 
the solar-system test, particularly the test in the atmosphere. 

• z ~ 0(1): The deviation from GR needs to be significant in order to explain the cosmic 

acceleration at the present epoch. 

According to the above constraints, in the viable f(R) models of the late-time cosmic 
acceleration the deviation from GR should be small when z > 70 but become significant at 
the recent epoch. The above constraints give simple, clear requirements one can conveniently 
utilize to examine the viability of the f(R) models with various functional forms of f(R). 
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